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I.  INTRODUCTION 


A.  Background 

There  Is  a  logical  sequence  of  events  In  the  design  of  any  defensive 
weapons  system.  Initially,  It  Is  recognized  that  a  known  or  potential  threat 
exists  for  which  one  or  more  countermeasures  are  desired.  There  are  specific 
threat  capabilities  that  existing  systems  cannot  defeat,  necessitating  either 
modification  of  the  existing  system  or  the  development  of  a  new  one.  Three 
groups  of  people  -  strategists,  designers,  and  analysts  -  participate  In  the 
development  and/or  modification  process. 

Initially,  strategists  and  planners  develop  the  general  concepts  of  a 
weapons  system  to  defeat  the  threat.  Their  decisions  are  based  on  knowledge 
of  flight  capabilities  of  the  threat  obtained  from  Intelligence  reports  and 
other  sources.  Defensive  system  specifications  are  established  to  ensure 
defeat  of  the  threat  while  minimizing  Its  Impact  on  military  operations. 
Designers  then  assemble  hardware  In  accordance  with  system  specifications. 
Analysts  have  a  two-fold  role  In  the  process:  (1)  they  develop  techniques  for 
testing  and  comparing  alternate  designs  and  (2)  they  translate  the  results  of 
simulation  studies  Into  usable  design  criteria  for  hardware  designers.  It  Is 
this  role  which  provided  the  Impetus  for  the  work  reported  here. 

This  report  addresses  one  particular  facet  of  an  analyst's  contribu¬ 
tion:  namely,  the  development  of  an  algorithm  for  calculating  and  plotting 
Inner  tracking  limits  of  a  tracking  system.  This  work  Is  an  extension  of  the 
work  of  Campbell  and  Christensen  (Report  No.  RL-TN-70-3,  "Tracking  Mount  - 
Relationships  for  a  Constant  Velocity  Maneuvering  Target").  The  notation  and 
basic /klnematlcal  relationship^  are  retained. 

The  term  "tracking"  refers  to  the  process  of  positioning  a  rigid  body 
so  that  for  all  Instants  of  time  a  line  in  the  rigid  body  passes  through  a 
moving  point.  Common  examples  of  such  bodies  are  guns,  guidance  seekers,  and 
missile  launchers.  The  rigid  body  Is  typically  constrained  to  rotate  about 
horizontal  and  vertical  axes  (elevation  and  traversing).  The  torques  and/or 
forces  necessary  for  tracking  are  Interrelated  functions  of  object  point 
movement,  rigid  body  mass  distribution,  rigid  body  kinematics,  and  friction. 
Since  It  Is  Impossible  to  design  a  system  to  track  for  all  positions  and 
velocities,  some  knowledge  of  the  limitations  and  capabilities  of  a  proposed 
tracking  system  Is  necessary  In  the  Initial  design  phase.  Equally  Important 
Is  the  ability  to  estimate  the  locus  of  points  In  space  where  a  tentative 
tracking  system  can  be  expected  to  track  a  specified  moving  point. 

Estimates  of  the  rotational  speeds  and  accelerations  essential  to  an 
acceptable  system  must  be  available  before  design  can  begin.  Obviously,  these 
data  are  Influenced  by  the  kinematics  of  the  moving  point  and  significantly 
Impact  the  choice  of  system  driver(s).  Equations  for  calculating  these  data 
as  functions  of  point  position  and  velocity  were  derived  as  part  of  this 
effort  and  can  be  used  Independently  of  the  plotting  segments  of  the  program 
to  check  the  adequacy  of  a  proposed  design. 
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The  referenced  report  of  Dean  and  Campbell  considered  the  same  problem 
treated  here,  but  their  procedure  yielded  larger  regions  where  tracking  Is 
possible  than  those  generated  when  more  realistic  characteristics  of  the 
driver  are  considered.  It  Is  more  realistic  and  more  accurate  to  examine  each 
point  In  space  In  relation  to  the  driver  angular  velocity  and  angular  accel¬ 
eration  associated  with  a  target  passing  through  the  point  to  see  If  the  com¬ 
bination  of  velocity  and  acceleration  are  within  the  capabilities  of  the 
system  driver.  Velocity  and/or  acceleration  bounds  may  preclude  tracking, 
l.e.,  tracking  may  require  an  angular  velocity  In  excess  of  maximum  driver 
speed  or  a  torque  larger  than  the  driver  Is  capable  of  supplying. 

This  report  addresses  this  essential  activity  which  resulted  In  an 
algorithm  for  establishing  the  Inner  tracking  boundaries  of  a  tracking  mount. 
The  algorithm  was  Implemented  on  an  ALPHA-MICRO  computer  system.  A  finite 
number  of  parallel  constant  speed  and  constant  altitude  paths  are  examined  to 
Identify  points  along  the  paths  where  necessary  angular  accelerations  and 
velocities  are  not  within  the  torque  and  speed  capabilities  of  the  driving 
mechanisms.  The  output  of  the  program  Is  a  plot  of  those  points  on  the  curve 
formed  by  the  Intersection  of  the  surface  bounding  points  In  space  where 
tracking  la  not  possible  (Appendix  B).  The  analysis  discussed  here  Is  limited 
to  an  elevating  mechanism.  Modification  of  the  program  to  Investigate  a  tra¬ 
versing  mechanism  would  be  straightforward. 

This  work  consisted  of  developing  computational  algorithms  and  digital 
computer  codes  for  calculating  and  plotting  boundaries  which  partition  the 
space  surrounding  a  cracking  mount  Into  disjoint  regions  where  cracking  Is 
possible  and  where  tracking  Is  Impossible.  The  limitations  are  based  on  an 
assumed  linear  relationship  between  the  angular  velocity  and  torque  capacity 
of  elevating  and  traversing  motors  (Appendix  A).  A  maximum  angular  velocity 
for  Che  driver  Is  Implied;  l.e.,  Che  magnitude  of  the  angular  speed  In  any 
direction  Is  bounded.  A  point  on  a  tracking  boundary  Is  determined  If  Che 
required  torque  exceeds  that  available  from  the  driving  motor  and  the  accom¬ 
panying  angular  rate  has  a  magnitude  and  direction  which  exceeds  the  maximum. 

The  output  of  the  program  Is  a  closed  curve  In  a  plane  of  constant 
elevation.  The  shape  of  Che  curve  Is  a  function  of  target  speed,  target 
elevation,  and  driver  characteristics.  Points  belonging  to  the  region 
enclosed  by  the  curve  represent  points  In  space  for  which,  under  the  con¬ 
ditions  specified,  a  target  cannot  be  tracked.  If  suitable  graphics  software 
Is  available,  several  plots  for  different  altitudes  can  be  plotted  on  a  single 
sheet  simulating  a  three-dimensional  plot  of  tracking  boundaries. 
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II.  MATHEMATICAL  BACKGROUND 


An  xyz-coordlnate  system  was  chosen  where  the  y-axls  Is  parallel  to 
the  flight  path  of  an  approaching  target,  with  the  positive  y-axls  directed 
toward  the  oncoming  target.  The  positive  z-axls  Is  upward,  and  the  positive 
x-axls  Is  chosen  to  form  a  right-handed  coordinate  system  (Fig.  1). 


Figure  1.  Illustration  of  an  xyz  -  coordinate  system. 
Definitions 


^  -  mount  to  target  heading  (azimuth) 

e  ■  mount  to  target  altitude  (elevation) 

a  *  target  x-coordlnate  (standoff  distance) 

c  -  target  altitude 

V  «  target  speed 

Using  Figure  1,  It  Is  easy  to  establish  the  geometric  relationships 

X  "  r  cos9  slni^ 
y  “  r  cos 9  cos(j, 

where  r  Is  the  distance  from  the  origin  to  the  moving  point.  Under  the 

conditions  of  the  problem  the  altitude,  c.  Is  given  by  the  equation 

c  ■  r  cos 9.  (1) 

From 

y  -  -  V  (2) 
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and 


x/y  ■  tan^ 

(3) 

It  follows  that 

•  2  2 
-  (xv/y  )cos  ♦ 

(4) 

2 

■  (v/a)sln 

(5) 

Also 

•  2 

9  «  (v/c)co8(^8in  9. 

(6) 

Rquaclona  (5)  and  (6)  can  be  differentiated  to  obtain  expressions  for  angular 
accelerations  as  functions  of  position  onl7. 

Consider  6.  It  Is  given  by 

e  ■  -  (iv/c)sln^sln29  +  2( 6v/c)coS(^slnecose  .  (7) 

Any  rigid  body  experiencing  an  angular  acceleration  must  be  acted  upon  by  a 
torque  about  the  axis  of  rotation  given  by  the  relationship 

T  -  10  (8) 


where  I  Is  the  moment  of  Inertia  of  the  mount  with  respect  to  the  axis  of 
rotation.  For  tracking,  this  torque  must  be  less  than  the  maximum  torque  out¬ 
put  available  for  the  associated  angular  velocity.  Driving  motor  charac¬ 
teristics  (Appendix  A)  dictate  that  the  maximum  torque  available  for  a  given 
angular  velocity  9  Is  given  by: 


18  -  T  -  Tg  -  k9.  (9) 

Substituting  the  expressions  for  9  and  9  from  equations  (6)  and  (7)  above  and 
performing  simple  algebraic  manipulations  yield  the  conditions 


2  2  2  2  2 
f(9,(>)  ■  (v  l/c  Tg)sln  9r-sln  (^tan9  +  2cos  (^cos9sln9 


+  (ck/vI)coS(>l  ■  +1 


as  necessary  to  establish  a  point  on  the  tracking  boundary.  If  |f(9,(t))  |<  1, 
the  point  In  the  plane  z  -  c  specified  by  the  three  parameters  9,  (ji,  and  c  is 
a  point  In  space  for  which  tracking  Is  possible.  Pairs  (9,()))  for  which 
I  f(0,^)  |>  1  locate  points  In  the  plane  z  ■  c  where  tracking  Is  not  possible. 
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Before  describing  the  procedure  to  establish  tracking  boundaries,  some  of 
the  Important  properties  of  the  function  f  are  required.  It  Is  instructive  to 
plot  the  function  f  for  fixed  (j*  and  for  9  varying  between  0  and  (Fig.  2). 
With  c  also  fixed,  this  variation  In  9  effectively  generates  points  along  the 
line  of  Intersection  of  the  planes  z  «  c  and  ■  constant,  and  the  ordered 
triple  (0>^,c)  locates  a  unique  point  along  this  line.  Observe  that  for  small 
there  are  two  values  of  9  satisfying  f(0,^)  ■  1  and  one  value  of  9  satisfy¬ 
ing  f(0,^)  ■  -1.  This  property  of  f  Introduces  a  complexity  Into  the  algo¬ 
rithm  for  determining  points  on  the  tracking  boundary  and  results  In  an  unex¬ 
pected  shape  for  the  area  defined  by  this  boundary. 


Consider  an  ordered  pair  (9,<^)  where  the  point  (9,f(9,i^))  belongs  to  the 
area  In  Figure  2  defined  by  the  Inequality 

I  fo.o;  I  <  1.  (11) 

Recall  that  the  ordered  pair  (9,(^)  defines  a  unique  point  (x,y,c)  In  the 
plane  z  •  c.  For  this  point  and  given  target  velocity,  9  and  9  can  be  calcu¬ 
lated  using  equations  (6)  and  (7)  above.  These  are  the  acceleration  and  velo¬ 
city  of  the  elevation  angle  necessary  for  tracking  a  target  passing  through 
this  point  In ^ space.  Then  a  point  In  Figure  A-1,  Appendix  A,  whose  coor¬ 
dinates  are  (9,19)  lies  between  lines  consisting  of  line  segments  AB  and  CD 
and  their  extensions.  If  In  addition 


then  the  given  driver  Is  capable  of  meeting  the  requirements  for  tracking  a 
point  moving  through  (x,y,c)  with  speed  v. 
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Figure  3  (valid  for  small  (^)  la  an  aid  In  Identifying  points  where 
tracking  Is  not  possible. 


Figure  3.  Tracklng/non-tracklng  In  plane  it  ■  constant. 

Line  segment  OD  denotes  the  intersection  of  planes  z  ■  c  and  it)  ■  Oq, 
where  i^q  Is  considered  small,  as  above.  The  three  angles  6^,  02,  9^^  are  the 
angles  defined  In  the  previous  paragraph.  Points  on  the  segment  RC  are  asso¬ 
ciated  with  values  of  6  satisfying 

<  9  <  02  .  (13) 

For  9  satisfying  this  condition,  f(9,^)  >  1  and  tracking  is  impossible. 
Also,  for 


9^  <  0  <  ir/2,  (14) 

<  '1  and  tracking  Is  Impossible  for  points  along  the  segment  OA.  The 
tracking  boundaries  for  fixed  altitude  consist  of  the  locus  of  points  A,  B, 
and  C  for  values  of  satisfying 

0  <  t  <  TT.  (15) 

Because  f(9,<^)  *  f(9t  ^  t)  the  tracking  boundary  Is  symmetric  to  the  plane 
X  *  0.  This  symmetry  Is  used  In  the  plotting  algorithm  described  below.  It 
can  be  shown  that  the  tracking  boundary  Is  not  symmetric  to  the  plane  y  >  0. 

The  algorithm  for  determining  boundaries  is  straightforward.  Parameters 
such  as  velocity,  height,  and  torque-angular  characteristics  are  assigned. 

The  values  of  the  function  f(9,4)  are  determined  for  different  values  of  9  and 
^  as  follows.  The  traversing  angle,  4,  in  Incremented  between  0  and  ir  radians 
starting  at  0  radians.  For  each  i^,  a  value  or  values  of  9  are  determined 
which  satisfy 

|f(0,o)|-+l.  (16) 
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These  angles  identify  limiting  torque-velocity  combinations  and  are  used 
to  define  a  point  on  the  tracking  boundary.  The  corresponding  rectangular 
coordinates  are  determined,  the  point  plotted,  and  the  process  repeated.  The 
procedure  is  straightforward  except  for  identifying  the  value  of  41,  say 
for  trhlch  two  values  of  9  exist  (Fig.  2).  This  value  of  41  is  the  upper  bound 
of  the  set  of  1^  for  which  there  are  two  ordered  pairs 

(9i.H).  1  -  1.2  (17) 

satisfying 

f(9,«)  -  1-  (18) 

For  ♦q  <  (^  <  IT  ,  there  are  no  real  values  of  9  which  satisfy 

f(9,*)  -  1,  (19) 

and  only  one  value  of  9  satisfying 

f(9,^)  -  -1.  (20) 

The  most  complex  part  of  the  computational  algorithm  is  the  logic  for 
recognizing  and  handling  the  special  situation  arising  when,  to  a  single  value 
of  (►,  there  are  two  values  of  9  for  which  f(9,(>)  ■  +1.  Ordered  pairs  (9,4>) 
satisfying  (19)  are  found  using  an  Increasing  sequence  of  values  of  9,4)  fixed. 
There  will  be  a  first  ^  for  which  there  is  no  9  satisfying  (19).  Let  this 
angle  be  called  (Fig.  2).  If  (>  is  decreased,  subsequent  values  of  9  in  the 
solution  pairs  (9,^)  of  (19)  will  increase  with  t/l  as  an  upper  bound.  The 
azimuth  angle  is  decremented  to  zero.  As  approaches  0  points  on  the 
tracking  boundary  approach,  the  z-axis  tangent  to  the  plane  is  x  -  0.  Finally, 
1^  is  Incremented  in  the  interval  with  the  search  for  9  satisfying  (20). 

It  should  be  clear  that  different  combinations  of  system  parameters  will 
result  in  different  tracking  boundaries.  In  addition,  many  combinations  of 
system  parameters  will  yield  the  same  plot.  For  example,  consider  the  factor 

v^I/c^T  (21) 

in  Equation  (10).  This  factor  is  dimensionless  and  every  choice  of  parameters 
yielding  the  same  value  for  this  ratio  will  result  in  the  same  tracking  boun¬ 
dary. 

The  intersection  of  the  tracking  boundary  and  a  plane  z  »  constant  is  a 
closed  curve  and  the  area  enclosed  will  approach  zero  with  Increasing  height. 
The  area  cannot  be  zero  for  a  finite  altitude  because  of  the  rapid  reversal  In 
9  necessary  for  small  standoff  distance  when  ^  is  approximately  90”.  As  a 
special  case,  if  the  standoff  distance,  a,  is  zero,  tracking  requires  an  Infi¬ 
nite  9  as  the  target  passes  directly  above  the  tracking  mount.  The  tme  de¬ 
rivative  of  the  elevation  angle  must  reverse  Instantaneously  to  track  a  point 
moving  directly  overhead. 
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Additional  insight  into  the  relationship  between  the  driver  requirements 
and  position  can  be  obtained  from  a  plot  of  a  torque-angular  velocity  rela¬ 
tionship  with  as  a  parameter  varying  between  0  and  t.  The  plot  will  begin 
at  the  origin,  since  no  torque  or  angular  velocity  is  required  for  points  at 
infinity.  The  larger  torque  and  angular  velocity  requirements  for  tracking 
close-in  targets  will  generate  points  more  distant  from  the  origin  and,  under 
some  conditions,  points  on  the  tracking  limit  boundary.  The  6  and  <{)  associ¬ 
ated  with  these  points  determine  a  point  on  the  spatial  tracking  boundary. 

The  boundary  surface  for  a  given  target  (v  •  constant)  can  be  determined  by 
considering  all  possible  standoff  distances  and  altitudes.  With  a  graphics 
system  capable  of  manipulating  and  displaying  3-D  geometric  data,  the  bounding 
surface  can  be  displayed  to  aid  in  visualizing  its  size  and  shape.  An  example 
obtained  using  an  Evans  and  Sutherland  PS-300  graphics  system  Is  shown  In 
Figure  4. 

As  discussed  above,  a  mount  will  lose  the  ability  to  track  when  the 
torque-angular  velocity  curve  crosses  the  tracking  limit  boundary.  Note  that 
there  are  two  limiting  phenomena,  those  imposed  by  the  sloping  straight  lines 
(torque)  and  those  imposed  by  the  vertical  lines  (velocity).  For  the  sloping 
lines,  the  Inability  to  track  is  characterized  by  a  deficiency  in  torque,  l.e., 
for  the  given  angular  velocity,  the  motor  Is  Incapable  of  generating  the 
required  torque.  The  vertical  lines  represent  upper  and  lower  bounds  on  the 
angular  velocities  of  the  motor.  The  driving  motor  is  Incapable  of  speeds 
outside  these  bounds. 

Selected  curves  to  Illustrate  several  possibilities  are  shown  in  Figure  5. 
The  ordinate  and  abscissa  are  dimensionless  variables  T/Tg  and  e/9nax’  Curve 
A  is  generated  by  a  target  flight  path  which  does  not  Intersect  the  Inner 
tracking  boundary  surface.  The  locus  of  points  representing  required  torque- 
angular  velocity  combinations  is  contained  within  the  closed  curve  specifying 
driver  characteristics.  There  are  many  sets  of  parameters  which  will  generate 
such  a  torque-angular  velocity  curve.  Large  values  of  "a"  and/or  "c"  or  small 
values  of  "I"  will  contribute  to  generating  similar  curves.  The  trajectory 
generating  curve  B  crosses  the  limiting  boundary  at  four  points.  Points  out¬ 
side  the  boundary  correspond  to  points  in  space  where  tracking  is  not  possible. 
It  can  be  concluded  from  curve  B  that  the  trajectory  which  generates  it  con¬ 
tains  two  intervals  where  the  target  is  Incapable  of  being  tracked.  The 
segments  1-2  and  3-4  lie  outside  the  characteristic  curve  and  correspond  to 
points  on  the  associated  trajectory  where  tracking  is  not  possible.  Curve  C 
is  generated  by  a  trajectory  for  which  velocity  as  well  as  acceleration  limi¬ 
tations  preclude  tracking. 


Figure  4. 


III.  SH0RTC(»1INGS  AND  DEFICIENCIES 


Throughout  this  report,  the  rotary  Inertias  have  been  assumed  Indepen¬ 
dent  of  position,  l.e.,  ’’I"  has  been  assumed  constant.  In  reality,  the 
moments  of  Inertia  are  time  and  position  dependent.  For  example,  the  moments 
of  Inertia  of  a  missile  launcher  are  a  function  of  the  number  of  missiles  In 
the  launcher,  and  the  moment  of  Inertia  of  a  mount  about  a  traverse  axis 
depends  on  the  elevation  angle  of  the  mount.  The  error  Introduced  by  assuming 
them  constant  will  depend  on  the  total  variation  of  I  over  the  elevation 
limits.  The  variation  of  I  with  elevation,  as  well  as  static  moments  of 
system  weights,  could  be  accounted  for.  Other  torque-speed  relationships  may 
be  more  realistic  than  the  linear  relationship  treated  herein. 

The  plots  Indicate  that  there  may  be  trajectory-elevation  combinations 
which  determine  disjoint  sets  of  points  where  tracking  Is  Impossible,  Indi¬ 
cating  that  the  ability  to  track  a  specific  target  can  be  lost,  regained,  and 
lost  again  before  being  acquired  for  the  last  time.  (If  the  range  Is  large 
enough,  any  system  Is  capable  of  tracking  a  moving  target,  assuming  adequate 
electronic  acquisition  and  tracking  capability.)  This  Information  Is  of 
little  value  unless  the  control  system  Is  capable,  after  losing  the  ability  to 
track,  of  moving  the  tracking  mount  to  the  next  available  point  to  resume 
tracking  when  the  target  reaches  the  point.  From  a  practical  point  of  view, 
the  most  useful  knowledge  available  from  the  algorithm  Is  the  location  of  the 
target  when  tracking  capability  Is  first  lost. 

IV.  CONCLUSIONS 

The  techniques  described  here  are  useful  for  visualizing  the  regions  In 
space  where  a  given  tracking  mount  can  track  a  target.  They  are  too  time  con¬ 
suming  to  be  used  In  a  real-time  situation  and  are  not  suitable  as  a  basis  for 
electing  alternate  procedures  when  tracking  becomes  Impossible.  However, 
through  the  use  of  predefined  tables,  the  location  of  the  next  point  at  which 
tracking  Is  possible  can  be  determined  If  tracking  capability  is  lost.  A 
control  algorithm  can  be  designed  to  position  the  mount  to  reacquire  the 
target  at  some  subsequent  position  when  tracking  again  becomes  possible. 

This  work  demonstrates  the  possibility  of  creating  a  plot  of  the  Inner 
tracking  boundaries  for  a  simple  tracking  mount  tracking  a  constant  velocity 
target.  The  techniques  can  be  extended  to  treat  more  complex  situations  where 
the  target  changes  velocity  or  rotary  Inertias  are  functions  of  position. 

The  acceleration  and  velocity  functions  can  be  used  during  preliminary 
design  for  sizing  motors  to  meet  minimal  ranges  of  engagement. 
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APPENDIX  A 

MOTOR  CHARACTERISTICS 


The  aflsuved  torque-speed  relationship  for  the  elevation  and  traversing 
drivers  Is  shown  In  Figure  A-1.  The  two  Inclined  lines  express  the  output 
torque  as  a  function  of  rotational  speed  for  a  constant  applied  external 
voltage*  Note  for  each  angular  speed  there  are  two  torques  possible  depending 
on  direction  of  rotation  and  the  sign  of  the  voltage.  One  corresponds  to  a 
situation  In  which  the  applied  voltage  creates  a  torque  In  the  direction  of 
rotation.  The  other  describes  the  situation  when  the  torque  acts  as  a  brake, 
providing  an  angular  acceleration  opposite  to  the  direction  of  rotation  and 
reducing  the  speed  of  rotation.  Line  AB  graphs  the  relationship  when  the 
external  voltage  Induces  rotation  In  one  direction.  Line  CD  graphs  the  rela¬ 
tionship  If  voltage  Is  reversed.  Therefore,  for  a  given  angular  speed  there 
can  be  two  output  torques,  depending  on  the  sign  of  the  applied  voltage. 


6 (ANGULAR 
VELOCITY) 


Figure  A-l.  Torque-speed  relationship  for  elevation  and  traversing  drivers. 

It  Is  further  assuaed  that  If  the  magnitude  of  the  applied  voltage  Is 
changed,  the  area  enclosed  by  the  characteristic  curves  changes  proportionate¬ 
ly.  For  example;  If  the  magnitude  of  the  applied  voltage  Is  halved,  the  maxi¬ 
mum  magnitude  of  the  angular  speed  Is  also  halved,  reducing  by  a  factor  of 
four  the  torque-angular  velocity  combinations  accessible  to  the  control  sys¬ 
tem.  The  coordinates  of  each  point  In  the  plane  are  a  pair  of  numbers  repre¬ 
senting  the  mechanical  state  of  a  motor.  Those  points  In  the  Interior  region 
represent  physically  possible  states,  i.e.,  through  the  use  of  appropriate 
control  algorithms  It  Is  possible  for  the  motor  to  operate  at  those  states. 

At  the  given  rotational  speed  It  Is  possible  to  adjust  the  applied  voltage  so 
that  the  required  torque  Is  available.  Combinations  of  torque  and  angular 
velocity  associated  with  points  In  the  exterior  are  physically  Impossible,  and 
target  velocities  and/or  positions  yielding  such  combinations  are  positions 
for  which  tracking  Is  Impossible. 

It  Is  assumed  that  any  combination  of  torque  and  speed  which  will  plot 
within  the  boundaries  of  the  parallelogram  can  be  met  by  the  elevating  or  tra¬ 
versing  driver;  I.e.,  a  control  system  can  be  designed  to  meet  the  torque- 
speed  requirement  represented  by  point  P.  The  figure  above  shows  that  a 
requirement  represented  by  P  cannot  be  met.  The  torque  associated  with  point 
P  is  larger  than  can  be  provided  by  the  driver  at  the  corresponding  angular 
velocity. 
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APPENDIX  B 


ALGORITHM  FOR  TRACKING  BOUNDARIES 
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The  algorithm  for  calculating  and  plotting  tracking 
boundaries  was  Implemented  on  an  ALPHA  MICRO  minicomputer 
Interfaced  with  a  Hewlett-Packard  digital  plotter.  The 
programming  language  la  BASIC,  but  the  program  can  be  converted 
with  minimal  difficulty  to  some  other  language. 


THIS  PROGRAM  DETERMINES  THE  TRACKING  LIMITS  IMPOSED  BY  LIMITATIONS 
OF  THE  DRIVERS  OF  ELEVATING  AND  TRAVERSING  MECHANISMS. 

IT  IS  ASSUMED  THAT  THERE  IS  A  LINEAR  RELATIONSHIP  BETWEEN  THE  TORQUE 
OUTPUT  AND  ANGULAR  VELOCITY  OF  A  DRIVING  MOTOR  OF  THE  FORM; 

T  -  (+/-)  TS  +  K  *  0,  WHERE 
T  -  TORQUE, 

TS  -  STALL  TORQUE, 

K  -  SLOPE  OF  THE  TORQUE-ANGULAR  VELOCITY  LINE, 

0  -  THE  ANGULAR  VELOCITY. 

THE  ABSOLUTE  VALUE  OF  THE  ANGULAR  VELOCITY  MUST  REMAIN  LESS  THAN  TS/R. 

REM  ******  DEFINE  PARAMETERS  AND  DIMENSIONLESS  VARIABLES  ****** 

SIGNIFICANCE  11  :  REM  ESTABLISHES  FIXED  POINT  ACCURACY 
GOSUB  600  :  REM  READ  DEFAULT  VARIABLES  AND  INTIALIZE  PLOTTER 
CALL  HP'INI 

REM  ****  MENU  **** 

MENU;  PRINT  TAB(-1,0) 

?“  MENU  OF  MTHETA  “ 

?;?;? 

?"  1-NEW  DATA" 

?"  2-NEW  PLOT" 

?"  3-RUN  ACCEL.  PLOT" 

?"  4-RUN  VELOCITY  PLOT" 

?"  5-END" 

?"  6-GRID" 

^  MtS 

INPUT"", NUM 
?TAB(-1,0) 

ON  NUM  GOTO  500,300,100,150,3000,900 
GO  TO  MENU 


REM  ****  PLOT  LIMITS  BASED  ON  TORQUE  **** 

100  DT  -  DTI ;B-1; GOSUB  200 
DT  -  DTI ;B«-1; GOSUB  200 
GOTO  MENU 

rEM  ****  PLOT  LIMITS  BASED  ON  VELOCITY  **** 


B-1 


150 


ZZ-1/DR1;IF  ZZ  >  1  THEN  ?  TAB(-1,0):POR  I-l  TO  10:?;NEXT  I;& 

?"N0  BOUNDARY  LIMITS  EXIST  FOR  THIS  CASE!!";FOR  I-l  TO  2000: NEXT  I:& 
GOTO  MENU 

T-ASN(SQR(ZZ)):  CP-1:SP-0:B-1: ?  PON$;CALL  PL’T:?  PD$ 

PH-0:GOSUB  400 

PH-PI-PH:CP-COS(PH):SP-SIN(PH):?PON$:CALL  PL’T:?  PD$ 

GOSUB  400 

T-ASN(SQR(ZZ));CP-1:SP-0;B— 1;?  PON$:CALL  PL’T:?  PD$ 

PH-0:GOSUB  400 
GOTO  MENU 

REH  ******  PLOTTING  ROUTINE  FOR  BOUNDARY  LIMITS  BASED  ON  ATTAINABLE  TORQUES 

REM  ******  INITIALIZE  AZIMUTH  AND  ELEVATION(PH  AND  T) 

EEM  ******  LOOK  POK  VALUE  OT  THETA  WHICH  MAKES  ABS(AC)-1 

200  T-0;PH-0 

GOSUB  1100:CP-1;SP-0:?PON$:?PU$:CALL  PL’T:?PD$ 

REM  ******  IF  AC>0  THERE  ARE  THREE  VALUES  OF  THETA  FOR  WHICH  ABS(AC))-1 
KEM  ******  THERE  ALSO  EXISTS  A  PH(PH-PH*)  FOR  WHICH  THERE  ARE  ONLY  TWO  ROOTS 
REM  ******  THIS  OCCURS  IF  THERE  EXIST  THETA  SUCH  THAT  AC  -  1. 

IF  AC  >  1  THEN  T-T-DTI;FLAG  -1:GOTO  1200  :  REM  CASE  FOR  THREE  ROOTS 

FLAG  -  0  :  REM  IDENTIFIES  CASE  HAVING  ONLY  ONE  ROOT  FOR  GIVEN  PHI 
PH  -  0  :  GOTO  1305 

EEM  ****  PLOT  VELOCITY  IMPOSED  BOUNDARIES  **** 


400  PH-PH+DP:CP-COS(PH);SP-SIN(PH) 

420  ZZ1-ABS(ZZ/CP);IF  ZZ1>1  THEN  ?  P0$; ?POFF$: RETURN 

440  T-ASN(SqR(ZZl)) 

CALL  PL’T:  GOTO  400 


1000  REM  ******  ROOT  FINDING  ALGORITHM  ****** 

1100  REM  *****  FIND  SMALLEST  THETA  FOR  WHICH  ABS(AC)-1 
1105  S-SIN(T);CT-COS(T) 

AC-DR3*S*S*(2 . *S*CT+0R2 ) 

IF  ABS(AC)  <  ITHEN  T^+0T:G0T0  1105 

IF  ABS(AC-SGN(AC))<.01  THEN  ?PON$;CALL  PL ’T:?PD$: RETURN 
T^-DT:DT-DT*.1;T-T-H)T:G0T0  1105 

1200  REM  *****  CASE  FOR  THREE  ROOTS  ****** 

REM  *****  ENTER  THIS  ROUTINE  WITH  PHI (AZIMUTH)  FIXED  AND  SMALLEST  THETA  FOR  AC-1 
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RBM  *****  GET  PRECISE  VALUE  OF  THETA  AND  PLOT  INITIAL  POINT  ****** 


t-t-dt 

1205  DT-DTI :  PH-PH-H)P 

CP-COS(PH):SP-SIN(PH) 

1210  S-SIN(T):CT-COS(T) 

AC-DR3*S*S*(-SP*SP*S/CT+2.*CP*CP*S*CT+DR2*CP) 

IF  AC  <0  THEN  GOTO  1248: REM  THETA  ROOT  DOES  NOT  EXIST 
IF  AC  <1THEN  T-T+DT'.GOTO  1210 

IF  (AC-1)<.01  THEN  TT-T:CALL  PL'T: ?PD$:GOTO  1205 
T-T-DT :DT-.l*DT:T-T+OT: GOTO  1210 

REM  DECREMENT  PH  UNTIL  PH<DP 
1248  T  -  TT 

1250  DT-DTI :PH-PH-DP: IF  (  PH  <  DP  )  THEN  GOTO  1300 
CP-COS (PH) : SP-SIN(PH) 

1255  S-SIN(T):CT-COS(T) 

AC-DR3*S*S  * ( -SP*S P*S /CT+2 . *CP*CP*S  *CT+DR2 *CP ) 

IF  AOl  THEN  T-T-H)T:IF  T  >  PI/2  THEN  GOTO  1300:  ELSE  GOTO  1255 
IF  1-AC  <  .01  THEN  CALL  PL’T -.GOTO  1250 
1260  T-T-DT :DT-.1*DT:T-T+DT: GOTO  1255 


1300  T-PI/2-DT 

1305  DT-DTI :PH-PH+DP: IF  PH>PI/2  THEN  GOTO  1350 

CP-COS(PH) : SP-SIN(PH) 

1310  S-SIN(T):CT-COS(T) 

AC-DR3*S*S*(-SP*SP*S/CT+2 .*CP*CP*S*CT+DR2*CP ) 

IF  AC  <  -1  THEN  T  -  T-DT:GOTO  1310 
IF  AC+1  <  .01  THEN  CALL  PL’T:GOTO  1305 
T-T+0T;DT-.1*0T: T-T-DT: GOTO  1310 

rEM  ****  LET  PH  INCREASE  TO  PI 
1350  T-T-DT 

1355  DT-DTI :PH-PH-K)P; IF  PH>PI-DP  THEN  GOTO  1400 
CP-COS(PH) : SP-SIN(PH) 

1360  S-SIN(T):CT-COS(T) 

AC-DR3*S  *S* ( -SP*SP*S / CT+2 . *CP*CP*S  *CT-H)R2  *CP ) 

IP  AC  >  -1  THEN  T-T+DT;IF  T  >  PI/2  THEN  GOTO  1380:  ELSE  GOTO  1360 
IF  -1-AC<.01  THEN  CALL  PL’T:GOTO  1355 

1380  T-T-DT ;DT-.1*DT: IF  DT  <  lE-10  THEN  GOTO  1400  ;  ELSE  T-r+DT:GOTO  1360 

rEM  ****  LET  PHI  DECREASE  TO  DP  LOOK  FOR  THETA  SUCH  THAT  AC  -  1 

1400  IP  FLAG  -  0  THEN  GOTO  1600:  REM  CASE  FOR  WHICH  AC  DOES  NOT  EQUAL  1 
DT-DTI :PH-PH-DP 
CP-COS  (PH)  :SP-SIN(PtI) 

1410  S-SIN(T):CT-COS(T) 

AC-DR3*S*S* ( -SP*SP*S /CT+2 . *CP*CP*S  *CT+DR2*CP ) 

IF  AC  <  1  THEN  T-T-DT: IF  T  >  0  THEN  GOTO  1410 

IF  T  <  0  THEN  GOTO  1600: REM  MO  VALUES  OF  THETA  EXIST 
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1>T-«)T;PH-PH+0P 
1430  DT-DTI : PH-PH-DP 

CP-COS(PH):SP-SIN(PH) 

1440  S-SIN(T):CT-C0S(T) 

AC-DR3*S*S*(-SP*SP*S/CT+2 . *CP*CP*S*CT+DR2*CP) 

IF  T  <  0  THEN  GOTO  1500: REM  NO  THETA  CAN  BE  FOUND 
IP  AC  <  1  THEN  T-  T-DT;GOTO  1440 
IP  AC  -  1  <  .01  THEN  TT-T:CALL  PL’T:GOTO  1430 
T»^+OT:DT-.1*DT:T-T-DT:GOTO  1440 

REM  ****  increase  phi  TO  PI 

1500  T-TT 

1530  DT-DTI :PH-PH-H)P: IP  PH>PI  THEN  GOTO  1600 

CP-COS(PH):SP-SIN(PH) 

1550  S-SIN(T):CT-COS(T) 

AC-DR3*S*S*(-SP*SP*S/CT+2 . *CP*CP*S*CT+DR2*CP ) 

IP  AC  >  1  THEN  T-T-DT:GOTO  1550 
IP  1-AC  <  .01  THEN  CALL  PL’T;GOTO  1530 
T-T*H)T:DT-.1*DT:T-T-DT:G0T0  1550 


REM  *****  PLOT  COMPLETED 

1600  ?PO$:?POPP$:REM  TURNS  PRINTER  OFF 
RETURN 

REM  THE  FOLLOWING  ARE  SUB  PROGRAMS  WHICH  CONTROL  THE  PLOTTER 
HP'PL: 

?  PRT 
PRT-"" 

RETURN 

HP 'INI: 

?  E$+".(" 

?  "IN;" 

?  "DP;" 

?  E$+".M:" 

?  E$+".H;;17:" 

?  E$+".N;19:" 

?  E$+".@:" 

?  E$+".)" 

RETURN 

300  INPOT"ORIGIN  COORDINATES  (X,Y)  IN  INCHES?  -  ";XXX,YYY:? 

330  INPUT"SCALE  FACTORS  FOR  X  AND  Y(UNITS  PER  INCH)?  -  ";SFXX,SPYY 

IP  SFYY  <-  0  CR  SPXX  <-  0  THEN  330 
? 

INPUT"  SELECT  PEN  0/1/2/3/4  - ";P 

? 

350  KR-1016 

YO-KK*YYY 
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XO-KK*XXX 

SPY-KK/SFYY 

SFX-KK/SFXX 


?  £$+".(" 

?  "SP"P";“ 

?  E$+".)" 

GOTO  MENU 

PL'T; 

R-C*COS(T)/SIN(T) 

YP-R*CP 

XP-B*R*SP 

XPI-INT(XP*SFX+'XO) 

YPI-INT ( YP*SFY+YO ) 

? "PA"STR(XPI) " , "STR(YPI ) “ ; " 
RETURN 

HP'END: 

?  E$+".(" 

?  "SP"0“;" 

?  E$+".)" 

RETURN 


REM  THIS  SUBROUTINE  DRAWS  THE  GRID 

900  ?TAB(-l,0) 

7-  *****  SELECTION  OF  GRID  PARAMETERS  *****":?:?:? 

?"  NOTE:  All  coordinates  of  the  gird  must  satisfy  0<"X<*15  and” 
?’'0<"Y<*11,  all  dimensions  Inches*”:?:? 

INPUT"ENTER  ORIGIN  IN  INCHES  -  X,Y  ”;X00,Y00 

INPUT"ENTER  SIZE  IN  INCHES  —  X,Y  ";PX,PY 

INPUT”ENTER  NO.  OF  DIVISIONS-  X,Y  ”;DPX,DPY 

INPUT”ENTER  PEN  NUMBER  0/1/2/3/4  ";P 

rEM  ******  CONVERT  TO  PLOTTER  COORDINATES  ****** 

XO-XOO*KK 

YO-YOO*KK 

DX-PX/DPX*KK 

DY-PY/DPY*KK 

X1-X0^•PX*^« 

Yl-YO+PY*KK 

?PON$ 

?  "SP"P";" 

rEM  ******  PLOT  GRID  ****** 

REM  ******  PLOT  HORIZONTAL  LINES  ****** 

Z-Xl-XO : YP-YO : N$»STR(DY ) 
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?PU$"PA"STR(XO)","STR(YO)";"  :  REM  MOVE  PEN  TO  ORIGIN 
920  ?PD$'*PR"STR(Z)",0;"PU$:Z— Z  :YP-YP+DT 

IP  Y1  >-  YP  THEN  ?"PRO,"N$";":GO  TO  920 

KEM  ******  PLOT  VERTICAL  LINES  ****** 


940  IF  Z  <  0  THEN  N$-STR(-DX) ;XP— XltXl— X0:& 
ELSE  N$-STR(DX);XP-XO 

Z-YO-Yl 

960  ?PD$"PRO,"STR(Z)";  "PU$;Z— Z;XP-XP+DX 

IP  XI  >-  XP  THEN  ?"PR"N$*,0;":  GOTO  960 

980  ?  E$+".)"  :  82M  TURN  OFF  PLOTTER 

GO  TO  MENU  :  REM  RETURN  TO  MENU 


REM  *****  MODIFICATION  OF  SYSTEM  DATA  ***** 

50Q  ?TAB(-l,0):?;?!? 

PRINT"INPUT  STANDOFF  DISTANCE  (";A;")  ";:INPUT"",  A 
PRINT-INPUT  TARGET  ALTITUDE  (";C;-)  -p.INPUT"",  C 
PRINT-INPUT  TARGET  VELOCITY  (-;V;-)  INPUT"-,  V 
PRINT-INPUT  MOTOR  STALL  TORQUE  (";ST;-) 

INPUT  --,ST 

PRINT- INPUT  SLOPE  OF  TORQUE  SPEED  CURVE 
INPUT  — ,K 

PRINT-INPUT  MOUNT  MOMENT  OF  INERTU  (-;MI;-) 

INPUT  —  ,MI 

PRINT-INPUT  INCREMENT  IN  THETA  FOR  PLOT  (";DTH;-) 

INPUT  — ,DTH:DTI-OTH/CF 

PRINT- INPUT  MAXIMUM  VALUE  FOR  THETA  (-;TTM;-) 

INPUT  ”,TTM:TMAX-TTM/CF 

PRINT- INPUT  VALUE  FOR  THE  ANGLE  PHI  (-;PHI;-) 

INPUT  — ,PHI:PHMIN-PHI/CF 

PRINT-INPUT  INCREMENT  IN  PHI  FOR  PLOT  (-;DPD;-) 

INPUT  --,DPD;DP-DPD/CF 
60SUB  650 
GOTO  MENU 

600  E$-CHR(27 ) ; PON$-E$+" . ( " : POFF$-E$+- . ) " ;PU$--PU ; " 

PD$--pD;-:CF-l80./3. 14159 

PI  -  3.141592654:K-1016:REM  CONV.  FROM  INCHES  TO  PLOTTER  COORDINATES 
READ  XXX,YYY,SFXX,SFYY,P,A,C,V,ST,K,MI,DTH,TTM,PHI,DPD 
DATA  1,1,1,1,0,1000,1000,1000,700,1140,450,5,90,0,1 

rEM  *****  CALCULATE  PARAMETERS  ***** 


650  DTI-DTH/CF : DP-DPD /CF : PHIMIN-PHI /CF ; TMAX-TTM/CF 

DR1-V*K/C/ST : DR2-K*C/V/MI : DR3-V*V*MI/C/C/ST 
RETURN 
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REM  ****  END  OF  PROGRAM  **** 
3000  CALL  HP 'END 
?TAB(-1,0) 

END 
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